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1. INTRODUCTION 

Switched reluctance machines (SRMs) are characterized by many advantages such as: a simple 
design, a wide range of possibilities for optimizing the operating modes when changing the speed and load (a 
possibility to adjust the output parameters of the electric drive), the overload torque capability, a high energy 
efficiency, reliability and low-cost, which are the reasons that make them considered as promising 
electromechanical devices [1], [2]. A complete summary of their performances, control methods and other 
capabilities can be found in [3]. The complexity of modern electric drive systems increases the need to use 
simulation methods for analysis and optimization based on various computer-aided design systems. Their use 
allows for determining the permissible deviations of the values of various parameters and control signals in 
electric drive systems from their calculated values so that the electric drive maintains its efficiency [4]. 

A powerful tool for an effective design, study and control of an SRM is modelling with a 
MATLAB/Simulink environment. SRM models are divided into two main categories-linear and nonlinear 
models [5]. Linear models are applicable when the motor assumes to remain magnetically unsaturated during 
operation. Furthermore, the phase inductances change linearly from the rotor position. Linear models of 
three-phase 6/4 SRM and three-phase 12/8 SRM are demonstrated in [5]-[7]. These models combine the 
ability to obtain acceptable results with minimal simulation time. However, the accurate analysis of SRM 
performance requires precision modelling of the behavior of this motor with its nonlinear magnetic 
characteristics. The built-in nonlinear models in the MATLAB/Simulink environment are types 6/4, 8/6 and 
10/8 SRM [8]. The nonlinear models of other types of SRMs are not included in the Simulink libraries. 

The three-phase 6/4 SRM has two magnetic poles per phase, while the three-phase 12/8 SRM has 
four. The electrical cycle of 6/4 SRM corresponds to a 90° rotation of the rotor position, while the electrical 
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cycle of 12/8 SRM corresponds to a 45° rotation. Besides, when the speed of 12/8 SRM is equal to half of 6/4 
SRM speed, the inductance cycles of both machines are equal. The turn-on and turn-off angles of the 12/8 
SRM are half of the same angles at 6/4 SRM. In addition, [9] discuss the possibility of modelling in Simulink 
the three-phase 12/8 SRM by modifying the existing model of the three-phase 6/4 SRM. 

Furthermore, the existing SRM model can be modified directly, or a new model can be added to the 
model library [9]. However, this approach does not allow to precisely model the behaviour of a specific type 
of 12/8 SRM. It should be noted that the flux linkage profile and the torque profile are changed with current 
and position. Besides, when operating at higher phase currents, the change with current is nonlinear. 

The method for structuring a nonlinear model of SRM considered in [10] allows the user to create 
their model with the ability to set multiple input parameters and implement nonlinear static characteristics. 
The initial data for the proposed model are the relationship of flux linkage as a function of position and the 
phase current. The necessary values of the current depending on the position and flux linkage and the torque 
depending on position and current are received, solving the corresponding differential equations. The use of 
partial derivatives should be avoided because they can become a source of error [5]. The waveforms obtained 
during simulations are not verified experimentally. 

Saidani and Ghariani [11], presented the results from the modelling of 12/8 SRM. The flux linkage 
curve is prepared through coupled to MATLAB FEMM simulation for different rotor positions and then is 
used to extract the other performance as static torque, inductance and co-energy. The look-up table, 
developed in the form of a block in the Simulink model, is obtained through MATLAB/Simulink. The 
simulation results are not confirmed with an experiment. The presented nonlinear models in [12]-[14] 
describe methods for structuring an SRM nonlinear model. These models use extensive analytical 
calculations and partial derivative calculations. Polynomial interpolation is also used. Thus, they emphasize 
the computational time of SRM simulations and the accuracy of the flux linkage partial derivative. These 
methods do not allow easy synthesis of an SRM model. 

The nonlinear model presented in [15] uses the data obtained from measurements of the static 
characteristics of the motor. They were implemented in the form of 2D tables in the MATLAB/Simulink 
environment. Disadvantages of the considered model are: a missing detailed description for structuring the 
model; a menu for setting the basic parameters of the studied motor; a method by which the data obtained 
from measurements are introduced into the model. The main disadvantages of the experimentally obtained 
static characteristics for magnetization and torque of the motor are the lack of precision in the measurement 
and the inability to stop the rotor at high currents without minimal slippage. To perform a simulation study to 
another type of SRM, it is necessary to provide the user with the capability to create their specific model 
based on catalogue data and the same time, save the computational time of SRM simulations. 

The present study presents a methodology for synthesising a custom nonlinear model of three-phase 
12/8 SRM in the Simulink environment. The procedure we used with MATLAB/Simulink is based on 
creating look-up tables that approximate the necessary relations. For this purpose, the approach given in [5] is 
applied. The required SRM characteristics of the tested 12/8 SRM are modelled by a FEM, as demonstrated 
in [16], providing us with the different magnetic data to realise the look-up tables. Thus, the used procedure 
allows avoiding all partial derivatives. The data for the look-up tables in excel files are entered using the 
source code for reading data in the Simulink software simulator. The individual blocks for structuring a 
proposed nonlinear model in MATLAB/Simulink are described step by step. 


2. MATHEMATICAL DESCRIPTION OF THE SRM NONLINEAR MODEL 
The expression of the phase voltage for a one phase of the SRM, describing the dynamic behaviour 
of the SRM, while assuming that there is no mutual coupling to other phases, is given by [1]: 


dipn(ipn Opn) 


Vpn (t) = ipn (t). Rpn z dt , (1) 


where Vpn is the voltage applied to the terminals of one phase, ipn is the phase current, Rpp is the resistance 
per phase, Apn (ipn Opn) is the flux linkage per phase at different excitation currents and rotor positions Oy). 
The flux linkage in an SRM phase is given by: 


Aon = Lpon(ipr Opn). iph (2) 


where Lyn (i, hw Onn) shows the dependence of the phase inductance on the rotor position and the excitation 
current. The torque produced by one phase Tpn can be approximated by: 


; iph 9A%pn(ipn 8 ph) 7. 
Ton (Con Opn) = P (3) 
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the total electromagnetic torque for the SRM with 3-phases is given by: 
Te = Xj-1 Toni (ipni Opri) (4) 


where j is the phase number. The following mechanical equations were used to obtain the output speed of an 
SRM: 


[B= Te -Troaa — Fo, (5) 
w = dð /dt, (6) 


where J is the machine inertia, Tioaq is the load torque, F is the friction coefficient and œ is the rotor speed of 
the SRM. 


2.1. The nonlinear SRM model 

This section presents the model synthesis using some of the techniques described in [10] and the 
embedded models in Simulink. The required static characteristics are obtained by the finite element method 
(FEM) [16]-[20]. The mathematical model of the SRM is divided into electrical and mechanical parts. The 
electrical part is represented by in (1) and (2), while the mechanical part is represented by in (3)-(6). 

Creating an adequate nonlinear mathematical model requires considering the highly nonlinear 
characteristics of the motor such as the inductance, the flux linkage and the torque. For this reason, these 
dependencies are presented in the two-dimensional tables (look-up tables) at different rotor positions and 
discrete values of the currents through the phases. The required SRM characteristics are built with the 
Infolytica MotorSolve software. The modern software for the finite element analysis FEM analysis is 
characterized with high accuracy and also includes an extensive database of materials from magnets, steels 
and conductors with measured values of properties provided by their manufacturers. The MotorSolve 
software simulates machine performance using equivalent circuit calculations and a unique automated finite 
element analysis machine [21]. Several analysis methods are included in this platform, which give the user 
control of the accuracy versus time. 


2.1.1. The SRM operational characteristics and inductance profile 

As the first step, a set of relations flux/current/position, current/torque/position, and inductance/ 
current/position are necessary to create the model. The graphical representation of the results for the flux 
linkages dependence Aan lions Onn) , obtained after the FEM analysis for the tested three-phase 12/8 SRM, is 
shown in Figure 1. The flux linkage values obtained by the FEM correspond to 90 rotor positions from 0° to 
45° (at resolution of 0.5°) and 100 different values of the current from 0 to 20 A (current resolution of 0.2 A). 

The next characteristic is the inductance profile Bei toni Op n)- Using FEM tools, the dependence of 
the inductance versus rotor position of the used three-phase 12/8 SRM for different excitation currents is 
derived. The obtained results for the phase inductance, displayed in Figure 2, correspond to 90 rotor positions 
from 0° to 45° (at resolution of 0.5°) and 99 different values of the current from 0.2 to 20 A (current 
resolution of 0.2 A). 

The most challenging task is the construction of the flux data Apn rem (ipw Onn) since the model 
calculates the current as a function of flux and position ipn (Apn_rem pn). Hence, it is necessary to process 
and rearrange the data Ayn rem (ipn Onn) to produce the required current typ (Apn_rem: Opn). The graphical 
representation for the set of curves ipn (Aon rem pn) in the Infolytica MotorSolve is presented in Figure 3. 
The obtained values correspond to 90 rotor positions from 0° to 45° (at resolution of 0.5°). It should be noted 
that the data for the flux linkages dependence App rem (ipn pn) and the inductance profile Lyn_rem (ipn Onn) 
corresponds to the values for one phase. Because the study assumes there is no mutual coupling between the 
phases, the look-up tables used to for one phase are also used to model all phases. 

The dependencies of the torque produced by one phase as a function of rotor position and phase 
current values Tyn (iw Opn) are shown in Figure 4. The obtained results correspond to 90 rotor positions 
from 0° to 45° (at resolution of 0.5°) and 100 different values of the current from 0 to 20 A (current resolution 
of 0.2 A). The obtained graphs and numerical results correspond to the operation of the SRM in motoring and 
generating mode. In addition, the results show the influence of the nonlinearity of the magnetization 
characteristic on the phase torque. 
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Figure 3. Infolytica MotorSolve 2D Figure 4. Infolytica MotorSolve 2D phase torque vs 
lyn (pn rem Onn) characteristics of the tested current characteristics of the tested three-phase 12/8 
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2.1.2. The implementation of the SRM nonlinear model in Simulink 

The schematic diagram of the SRM nonlinear model implemented in MATLAB/Simulink 
environment is presented in Figure 5. The block named source includes the conventional blocks used for the 
implementation of in (1)-(6) and the nonlinear characteristics of the motor, receiving from FEM. The 
simulation of each motor phase is realized using the blocks named phase (phase A, phase B and phase C). 
The approach in [10] is used. 


9 aseYd 


aounog 
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Figure 5. The structure of the nonlinear three-phase 12/8 SRM model 


Figure 6 shows the internal structure of the block named source. It contains the three lookup-tables 
denoted as 2-D lookup current (ipn (Apn_rem 8pn)), 2-D lookup inductance (Lon rem (ipw 8pn)), 2-D lookup 
torque (Typ (ipn 8pn)), the block for calculating the total torque and the rotor position block (position_rotor). 
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Figure 6. The block source schematic diagram 


The data for the two-dimensional tables (2-D lookup tables) in excel files are entered using the 
source code for reading data in the simulink software simulator in the PostLoadFcn menu. Snapshots of the 
PostLoadFcn menu with part of the code and the user interface, where the values for the current are entered, 
for the creation of the respective block 2D-Lookup table are shown in Figure 7 and Figure 8. 

The block denoted as Torque implements the mechanical in (4)-(6). Figure 9 illustrates their internal 
structures. The input data are the signals of the load torque TL and the total torque Te, which determine the 
rotor speed œ and the magnitude of the rotor position Teta. The internal structure of the rotor position block 
is shown in Figure 10. In this block, the rotating speed is integrated, and as a result, the angular position of 
the rotor for three phases is obtained. The block named phase A (phase B, phase C) includes a current source 
and a resistor connected in parallel, modelling the phase active resistance. The phase current, derived from 
the differential equation, is fed to input I of block phase A. Its output signal produces the voltage applied to 
the phase. 
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breakpoints3 = data1(2:end,1)'; Breakpoints 2: breakpoints2 
% Column indices for lookup table 
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% Output values for lookup table 
table_datal = data1(2:end,2:end); 
% Import the data from Excel for a lookup table 
data2 = xlsread('Torque’,'Sheet2'); 
% Row indices for lookup table 
breakpoints5 = data2(2:end,1)'; 
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Cancel Help Apply 9 [ox] | cancer Help Apply 
Figure 7. PostLoadFcn menu with source code part Figure 8. The user interface for entering FEM 
for reading data received from FEM analysis simulation results on the corresponding table 


The tested SRM is fed from the asymmetric bridge converter. This topology is the most widely used 
in SRM drives in terms of converter losses and cost, control performance, the independent phase windings, 
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allowing a fault-tolerant operation capability [7]. Hence, for example, the terminals named "A+" and "A-" 
(block phase A) connect the motor phase A to the corresponding bridge of the asymmetric converter circuit. 

The final synthesis of the three-phase 12/8 SRM nonlinear model includes the creation of a 
subsystem (block subsystem), shown in Figure 11, where the user interface (mask) sets the necessary input 
parameters such as a phase resistance (Q), a friction coefficient (N.m.s) and inertial value (kg.m.m). In 
addition, by checking the plot magnetization curves field of the subsystem, the flux linkages dependencies 
are displayed as shown in Figure 12. 

Figure 13 presents the schematic diagram used to simulate in MATLAB/Simulink the proposed 
nonlinear three-phase 12/8 SRM model. It consists of the following main blocks: block SRM 12/8 (nonlinear 
SRM model), block converter, block position_sensor, monitoring the position of the rotor and switching the phase 
transistors. The internal circuit of the phase bridge A of the asymmetric converter is shown in Figure 14. 
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Figure 13. The MATLAB/Simulink schematic diagram for the nonlinear three-phase 12/8 SRM model 
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Figure 14. H-bridge asymmetric converter — bridge phase A 


The input parameters are as follows: input voltage for motoring V=45 V and generating mode V=8 V; 
phase resistance Rph=1.795 Q, inertial value J=0.00057 [kg.m.m], friction coefficient F=0.0013 [N.m.s], the 
used MOSFETs are IRFP460 with the following parameters: drain-source on-resistance Rps(on=0.24mQ, 
body diode on-resistance Rd=0.01 ©, body diode voltage Vsp=1.8 V; the asymmetric bridge diodes are 
FR604-max. forward voltage Vr=1.3 V. Table 1 shows the specification of the tested three-phase switched 
reluctance motor type 12/8 model - HSSPWBKM- 1844. 


Table 1. Specification of tested three-phase 12/8 SRM 
Parameters of switched reluctance motor model HSSPWBKM-1844 


Parameters Symbol Value Units Parameters Symbol Value Units 

Number of phases Nph 3 - Inductance aligned position La 42.55 mH 
(measured at 0.2 A) 

Number of stator poles Ns 12 - Inductance unaligned Lu 5.22 mH 
position (measured at 0.2 A) 

Number of rotor poles Nr 8 - inertia J 0.00057 kg.m.m 

Phase resistance Rph 1,795 ohms friction coefficient F 0.0013 N.m.s 

Max. current per phase Iphmax 5:5: amperes Rated speed n 1500 rpm/min 

Phase voltage Vph 120 volts rotor arc angle Br 15.5 deg. 

Power delivered P 660 watts stator arc angle Bs 15 deg. 


3. RESULTS AND DISCUSSION 

In this section we proposed the performed simulation and experimental results. To verify the built 
nonlinear model accuracy, the series of simulation and experimental results are accomplished. The tested 
SRM is fed from the asymmetric bridge converter. Primarily, the results of the FEM analysis for the 
inductance profile are compared with those obtained experimentally. For this purpose, the current saturation 
method for SRM inductance measurement proposed in [22] is applied. The method uses low frequency and 
low voltage to minimize iron losses; no additional specific devices are needed; the setup is not complicated 
and compared to classic methods [23]-[27], a higher accuracy of the inductance profile measurement is 
achieved. Figure 15 shows the simulation and experimental verification of the inductance profile at various 
excitation currents for the tested three-phase 12/8 SRM. 
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Figure 15. Simulation and experimental verification of the inductance profile Lpp at various excitation 
currents for the tested three-pahse 12/8 SRM (HSSPWBKM- 1844) 


A very good alignment between simulation and experimental results is observed. In the present 
study, to confirm the accuracy of the built three-phase 12/8 SRM nonlinear model, the machine is energized 
using a voltage source control method. The operation of the motor is examined in quadrants I and II at motor 
speed n=1500 rpm. At the moment when the positive voltage (V+) feeds to one phase and the motor starts to 
rotate, the two MOSFETs (Q1, Q2) shown in Figure 14 must be on. Conversely, when the applied voltage 
has a negative value (V-), the two diodes (D1, D2) are used to ensure a continuous current mode [1]. The 
proposed nonlinear model is verified experimentally. Figure 16(a) presents a built three-phase asymmetric 
bridge converter with the controller. The measuring stand for experimental study of the operation of the 
three-phase 12/8 SRM with position and torque sensors is shown in Figure 16(b). 


Figure 16. Experimental study: (a) three-phase asymmetric bridge converter with the controller and 
(b) the measuring stand with position and torque sensors 


Depending on the operating mode and the direction of the motor rotation, there are different turn-on 
and turn-off angles of the motor phases. The oscillograms shown in the following figures correspond to 
motoring and generating mode, respectively, i.e. the operation in the first and second quadrants. Hence, the 
turn-on and turn-off angles must be selected in the rising and falling phase inductance for motoring and 
generating mode. Figure 2 shows that these areas are between 22.5° and 45° for first quadrant and 0° and 
22.5° for second quadrant. Therefore, to ensure a rapid increase and decrease of the phase current, it is 
necessary switching the applied phase voltage in these areas, leading angles corresponding to the minimum 
and maximum phase inductance (or 22.5° and 45°) for motoring mode and maximum and minimum phase 
inductance (or 0° and 22.5°) for generating mode. After comparing the voltages and currents of the tested 
SRM during operation in the I quadrant (motoring mode), the operation of the SRM in the II quadrant 
(generating mode) was checked by using a collector machine for the motor. During the simulation, load 
torque TL=0.65 Nm and the motor speed n=1500 rpm were specified. 
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Figure 17(a) illustrates from top to bottom the simulation waveforms of the one-phase voltage and 
the one-phase current at Oon = 25.3°, Ooff = 43.44°, torque load TL=0 Nm and speed n=1500 rpm for the 
motoring mode (I quadrant). Figure 17(b) illustrates the simulation waveforms of the same quantities - from 
top to bottom waveforms of the one-phase voltage and the one-phase current at Oon = 3.8°, Oo¢¢ = 20.944, 
torque load TL=0.65 Nm and speed n=1500 rpm for the generating mode (II quadrant). 


<Phase Voltage (V)> <Phase Voltage (V)> 
T T T T 
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T T T 


i i i i i 
0.52 0.525 0.53 0.535 0.54 0.245 0.25 0.255 0.26 0.265 


(a) (b) 


Figure 17. Simulation waveforms of the phase voltage and the phase current: (a) motoring mode at the motor 
speed n=1500 rpm, TL=0 Nm; and (b) generating mode at the n=1500 rpm, TL=0.65 Nm 


The following Figure 18(a) and Figure 18(b) present experimental waveforms of the one-phase 
voltage and the one-phase current in motoring mode and generating mode. Figure 18(a) illustrates for the 
motoring mode (I quadrant) from top to bottom the experimental waveforms of the one-phase voltage and the 
one-phase current at Oon = 24.13°, Oo¢¢ = 42.274°, torque load TL=0Nm and speed n=1500 rpm. 
Figure 18(b) illustrates the experimental waveforms of the same quantities for the generating mode (II 
quadrant)-from top to bottom waveforms of the one-phase voltage and the one-phase current at 
Oon = 3.63°, Oo¢¢ = 20.87°, torque load TL=0.65 Nm and speed n=1500 rpm. 

The next Figure 19(a) and Figure 19(b) present simulation waveforms for the motor speed at the 
operation in the first and second quadrants respectively. Figure 20(a) and Figure 20(b) illustrate the simulation 
results of the total torque Te at the operation of the SRM in motoring and generating mode respectively. The 
following Figure 21(a) and Figure 21(b) show the simulation results of the one-phase inductance at the 
operation of the SRM in the first and second quadrants respectively. Table 2 and Table 3 summarize the 
comparison between simulation and experimental results for the two operation modes. The relative error 
between simulation and experimental results is within 5-10%. The differences between them are mainly due 
to the inaccuracy in the measurement of the angular position, emerging electromagnetic processes, noise and 
vibration when the phase winding is excited. 


(a) (b) 


Figure 18. Experimental waveforms of the phase voltage and the phase current, (a) motoring mode at the motor 
speed n=1500 rpm, the torque load TL=0 Nm and (b) generating mode at the n=1500 rpm, TL=0.65 Nm 
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Figure 19. Simulation results of the motor speed at operation in: (a) the first quadrant and the torque load is 
TL=0 Nm and (b) the second quadrant and the TL=0.65 Nm 
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Figure 20. Simulation results for the total torque Te at operation in: (a) motoring mode at the motor speed 
n=1500 rpm, the torque load TL=0 Nm and (b) generating mode at the n=1500 rpm, TL=0.65 Nm 
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Figure 21. Simulation results of one-phase inductance: (a) first quadrant at the motor speed n=1500 rpm, the 
torque load TL=0 Nm and (b) second quadrant at the n=1500 rpm, TL=0.65 Nm 


Table 2. A comparison between simulation and experimental results in the motoring mode (I quadrant) 


Turn-on Turn-off Speed, Phase Phase Total torque 

angel ĝon angle Oop rpm voltage [V] current [A] Te[Nm] 
Experimental results 24.13 42.274 1500 48 2.05 0.31 
Nonlinear model simulation results 25.3 43.44 1500 45 2.02 0.305 
Relative error, % -4.85 -2.75 0 6.25 0.97561 1.61 


Table 3. A comparison between simulation and experimental results in the generating mode (II quadrant) 


Turn-on angel Turn-off Speed, Phase Phase Total torque 
Oon angle 6f rpm voltage [V] current [A] Te[Nm] 
Experimental results 3.63 20.87 1500 10 2 -0.059 
Nonlinear model simulation results 3.8 20.944 1500 9 22 -0.065 
Relative error, % -4.68 -0.35 0 10 -10 -9.83 
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4. CONCLUSION 

The paper describes a nonlinear three-phase 12/8 switched reluctance machine model suitable for 
use in a MATLAB/Simulink environment. To build the model, two-dimensional tables (look-up tables) 
obtained using the finite element analysis method and the Infolytica Motorsolve software simulator are used. 
The data for the look-up tables in excel files are entered using the source code for reading data in the 
Simulink software simulator. In order to ensure the model's reliability, 2D graphs of the simulation results 
from the FEM analysis and those generated from the MATLAB/Simulink are compared. A subsystem (block 
subsystem), where the user interface (mask) sets the necessary input parameters such as a phase resistance 
(Q), a friction coefficient (N.m.s) and inertial value (kg.m.m) was created. Furthermore, an additional block 
has been added in the proposed model to monitor the inductance profile during the simulation study. A series 
of simulation and experimental results are derived to verify the accuracy of the built three-phase, 12/8 SRM 
(HSS5PWBKM-1844) nonlinear model. A very good parity between simulation and experimental results is 
observed. The relative error between simulation and experimental results is within 5-10%. 
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